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Abstract 



A new analysis of the Lunar Prospector epithermal neutron data is presented, 
providing an improved map of the distribution of hydrogen near to the lunar 
poles. This is achieved using a specially developed pixon image reconstruction 
algorithm to deconvolve the instrumental response of the Lunar Prospector's 
neutron spectrometer from the observed data, while simultaneously suppress- 
ing the statistical noise. The results show that these data alone require the hy- 
drogen to be concentrated into the cold traps at up to 1 wt % water-equivalent 
hydrogen. This combination of localisation and high concentration suggests 
that the hydrogen is present either in the form of a volatile compound or as 
solar wind protons implanted into small regolith grains. 



1 Introduction 



Is there ice on the Moon? Finding the answer is important both for our under- 
standing of the Solar System and for planning future visits. The fact that the 
Moon's axis of rotation is inclined at only 1.6° from the normal to the ecliptic 
means that permanently shaded regions exist inside many craters within 10° 
of the lunar poles. These regions should be sufficiently cold ( < 100K) to pre- 



vent significant sublimation of water ice over billions of years (jVasavada et al. 



19991 1 . The flux of water reaching the Moon via impacts with comets and mete- 



orites and subsequently able to migrate to the relative safety of the polar 'cold 
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traps' should be sig n ifican t over the age of the solar system (jArnoldl . Il979 



Crider and Vondrakl . 120031 ; iButlerl . 119971 ). Hydrogen delivered by the solar 



wind can also create water by chemical re duction of the l unar regolith, provid- 
ing a nother potential source of polar ice (jArnoldl . 1 19791 ; I Crider and Vondrakl . 
2000h . 



Whil e the existence of ice on the Moon was considered in the 1960s (jWatson et al 
19611 ) . interest was rekindled in the 1990s, when two different methods were 



employed to address the issue, namely radar reflection experiments and neu- 
tron spectroscopy from Lunar Prospector. The radar method relies upon de- 
tecting the preferential reflection of the same sense of circular polarisation 
carried by the incident radiation. This coherent backscatter opposition effect 
could result from volu me scattering within a largely transparent medium, such 
as reasonably pure ice (iHapkd . Il990l ). However, lunar observations have shown 
that same-sense circularl y polarised radio waves are preferentially returned 



even from sunlit locations (IStacv et al.l . Il997l ; ISimpson and Tylerl . Il999l ; ICampbell et al. 



20061 ; ICampbell and Campbelll . 120061 ). making the most plausible explanation 



that this is surface scattering from blocky ejecta within craters, rather than 
volume scattering from large chunks of ice. 

The strongest evidence for water ice at the poles is provided by the deficit in 
the epithermal (0.3eV-100keV) neutron counts detected by the Lunar Prospec- 
tor. Cosmic rays impacting upon nuclei in the regolith liberate neutrons. These 
elastically and inelastically scatter off nuclei in the top few metres of the lunar 
regolith and some leak out of the surface before they suffer nuclear capture 
or spontaneously decay. Adding hydrogen to the regolith means that kinetic 
energy is more rapidly drained from the neutrons during elastic collisions. 
Consequently, fewer epithermal neutrons escape to be detected by orbiting 
neutron spectrometers. Thus, a deficit of epithermal neutrons is produced 
when there is a concen tration of hydrogen in the top ~ 1 metre of the regolith 



flFeldman et all Il99lh 



Lunar Prospector's neutron spectrometer detected a significant dip in the 
epithermal neutron coun t rate over both luna r poles, sugg esting an excess of 
hydrogen in these areas (IFeldman et al.l . Il998l . |2000| . 1200 ll ). These tantalising 
results did not definitively determine whether or not the hydrogen was confined 
to the cold traps. If the excess hydrogen were distributed throughout the 
polar regions, then it could just be solar wind-implanted protons bound in 
OH groups that were less likely to diffuse out of the locally cold regolith 
( Starukhina and Shkuratov . 2000l ). Even if the hydrogen were confined to the 
cold traps, then the temperatur es should be low enou gh that regolith grains 
will retain some free H atoms. IFeldman et al.l (120011 ) studied this in detail, 
concluding that the diffusivity of hydrogen was unlikely to be relevant within 
cold traps. When trying to discriminate between these various possible forms 
for the excess polar hydrogen, it would clearly be helpful to know whether or 
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not it is confined to the permanently shaded regions, and how concentrated 
this hydrogen is. 



2 Pixon Image Reconstruction 



Earlier studies of the Lunar Prospector data have used either the high-altitude 



100 km) epithermal data (IFeldman et all [1998) or sm oothed versions of the 



low-altitude (30 km) data (IFeldman et all |2000| . l200ll ). The response of the 
Lunar Prospector has a FWHM of 130 km and 45 km for high and low altitude 
orbits respectively. To determine whether or not the hydrogen is confined to 
the cold traps requires an approach that deconvolves the Lunar Prospector 
data while not introducing spurious complexity. This is the usual challenge in 
image reconstruction. Namely, given an input two-dimensional map, /, which 
is measured with an instrument that blurs it through a response function, B, 
and introduces some random noise, N, then the measured data satisfies 



D = I*B + N. 



The inclusion of noise means that the best that can be achieved is that the 
inferred image, /, avoids unnecessary complexity and produces a residual field, 

R = D — I * B, (2) 



that is statistically equivalent to the anticipated noise, N. 

The standard way to approach this problem is to consider the following equa- 
tion concerning conditional probabilities: 

p( iMD) = mlM) f\ M)P{M) - (3) 

p(D) 



M represents all aspects of the model for transforming the inferred image, J, to 
the data, D. This includes the assumed response function, B, the pixellisation 
of the data and the assumed characteristics of the noise. p(D) is a constant, 
and, to avoid introducing prejudice, so is p(M). Thus, 

p(I, M\D) oc p(D\I, M)p(i\M). (4) 



The term on the left-hand side of this proportionality is commonly known as 
the posteriori probability, P pos t = p{I,M\D), which should be maximised by 
the reconstruction. The two terms on the right are the likelihood (p(D\I, M) = 
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e -o.5x f or fae usua ] definition of the \ 2 goodness-of-fit statistic) and the image 
prior, Pip = p(I\M), respectively. The likelihood increases as the complexity 
of the reconstructed image, /, is increased and the fit to the data improves. 
In contrast, the image prior decreases with increasing complexity in I. This 
means that the optimal reconstruction (maximum P pos t) occurs at the point 
where increasing complexity in the reconstruction no longer improves the fit 
sufficiently to be justified by the information content in the data. 



A particularly suitable method for max imising the poster i ori probability is the 
pixon image reconstruction algorithm (jPina and Puetterl . Il993l : lEkd . l200ll ). A 
pixon can be thought of as a collection of pixels carrying a particular amount 
of information. Thus, the reconstruction of the Lunar Prospector data, which 
is well-sampled at the lunar poles and less so further away due to the polar 
orbit of the spacecraft, will have many small pixons to describe I near to 
the pole, and a few larger pixons where the signal-to-noise ratio decreases at 
lower latitudes. / is constructed by convolving a pseudoimage, defined in the 
pixel array and constrained to lie within physically acceptable count rates, 
with a position-dependent kernel, whose width represents the local pixon size. 
The optimum reconstruction is the one that combines simplicity of the recon- 
structed image with an acceptable residual field. 



2.1 Choosing the optimum number of pixons 



In order to determine the optimum number of pixons, n p , required to recon- 
struct an image, the total information in the data set is first determined via 



n n ix) 



<Wn n (x) 



(5) 



where n n (x) is the number of neutrons in pixel x, and v (x) = yn n (x)/a is the 
signal-to-noise ratio in pixel x. Note that the choice of a is somewhat arbitrary. 
How many a should a signal be before it counts as a piece of information? A 
value of a = 2 has been chosen, although no significant change is made to the 
reconstruction by choosing a — 1. 

To maximise the entropy of the reconstruction, the information placed into 
each pixon should be the same. Thus, for a given total number of pixons in the 
reconstruction, each should contain T/n p pieces of information. The Gaussian 
pixon, K^, was selected for each pixel such that 

v * K x = T/np. (6) 
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With the pixon distribution set for a particular total number of pixons, the 
algorithm then uses the Fletcher-Reeves conjugate gradient algorithm to min- 
i mise E r , the autocorrel ation of the reduced residuals (r = R/a) statistic 



(IPina and Puetterl . Il992l ). applied for lags of 1 pixel. Having thus found the 
best fit for this number of pixons, the posteriori probability, P pos t, can be cal- 
culated by multiplying together the likelihood (e _0 5x ) and the image prior, 
which is 

T! 



ip 





.(**p)\ 


rip 









(7) 



Using Stirling's approximation gives 

ln(P post ) w 0.5[lnT + (1 - n p )ln(27r) - n p m( T Ab) - X \ (8) 

Thus, the number of pixons that maximises ln(P post ) can be found, and this 
provides the optimal reconstruction. 

2.2 Coupled and decoupled reconstructions 



The reconstructed image, /, is defined by gathering signal from a pseudoimage, 
W . This gathering amounts to a convolution of the pseudoimage with the 
normalised pixon kernel appropriate for each pixel, namely 

/(x) = W * K x . (9) 



Constraints are placed on W, so that it cannot take unphysical values, such as 
count rates less than zero or greater than 22 epithermal neutrons per second, 
which is lar ger than the rate expe cted from regolith containing no hydrogen 



whatsoever ([Lawrence et all 120061 ). A set of 11 Gaussian pixons, with widths 
varying from 10 km to 250 km are available for the algorithm to choose. When a 
single pseudoimage is employed, the sunlit and permanently shaded pixels are 
treated equally. This is termed a coupled reconstruction. However, one should 
include some prior knowledge that sunlit pixels would not be expected to have 
a hydrogen content exceeding a few times that in the returned lunar samples. 
A loose limit of 15 counts per second is chosen for the minimum allowable 
epithermal rate in sunlit pixels. Introducing an independently determined map 
of permanently shaded pixels enables a decoupled treatment of the two types 
of pixel; sunlit and shadow. Each pixel type has its own pseudoimage, with 
different limits on their allowed count rates: sunlit - [15,22] and shadow - [0,22]. 
These pseudoimages only have non-zero elements for pixels of the type being 
considered. Thus, to avoid having the pseudoimage grow to compensate for 



5 



this near to the boundaries between pixel types, a mask, m type (x), of ones and 
zeros is set up for each pixel type, with the ones lying only in pixels of the 
chosen type. This mask is then convolved with the appropriate pixon kernel 
for each pixel to give a map of weights 

Mtype = Ultype * K x . (10) 



The resulting map, M type , has values lying in the range (0, 1] and can be used 
to define the decoupled reconstruction via 

f _ gW * gx) _ (n) 



In practice because the required pixon size is a continuous variable, whereas 
only a finite set of distinct pixon sizes are actually used, the reconstruction is 
created by interpolating between the nearest two available pixon kernels. By 
decoupling the shadow and sunlit pixel values, the effective number of pixons 
being used increases. To account for this, n p is replaced in the image prior 
calculation by n p G, where 

^pixel ^ ^typc( x ) 



with the sum being over all n P i xe i pixels x. 

Note that, for the decoupled reconstructions of the real data, the sunlit pixels 
take values in the ranges [18.5,20.5] counts per second for the north pole, and 
[18.3,20.7] counts per second for the south pole. For the north pole, this range 
corresponds to plausible amounts of hydrogen for occasionally sunlit terrain. 
The cons truction of the nort h and south pole shadow maps is described in 



detail by lElphic et al.l (120071 ). In brief, radar tomographic maps (Margot et 



al., 1999; U.S. Geol ogical Survey, 2002) were combined with the model of polar 



craters produced by lBussey et al.l (120031 ). Gaps were filled by interpolation and 



the whole resulting landscape was projected on a sphere of 1738 km radius, and 
illuminated with summer sunlight. Note that the permanently shaded areas 
are better determined around the north pole because the south pole shadow 
map was largely de rived from data taken during winter lighting conditions 



(lElphic et al.l . 120071 ) . The less well determined shadow map in the south may 
be behind the larger range of values being used by the sunlit pixels in this 
case. 

The left hand panels of Fig.l show the north (top) and south (bottom) pole 
epithermal data. These can be contrasted with the central and right hand 
columns in Fig.l, showing / from two different reconstructions of both data 
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sets. These four reconstructions all provide residual fields with acceptable \ 2 
values, and differ only because an additional shadow map prior was included 
for the decoupled reconstructions in the right hand column, unlike the coupled 
reconstructions in the central column. The coupled reconstructions are very 
similar to the count rate maps created by merely smoothing th e observational 



data u sing a Gaussian with FWHM of 40 km, as chosen by iFeldman et al 



( 120011 ) . The decoupling significantly changes the inferred hydrogen abundance 
within the cold traps, which are clearly visible because of the rapid count rate 
jumps across their edges. If a good fit is defined merely using x 2 , then there 
appears to be no need to impose the shadow map. The coupled reconstruction 
is simpler too, so the use of the shadow map can only be justified if the data 
demand it. 



3 Results 



The coupled and decoupled reconstructions differ significantly only in the ep- 
ithermal neutron count rates present in the cold traps. Thus, the residuals in 
the vicinity of these permanently shaded areas offer the best route to determine 
if the data contain sufficient information to discriminate between these two 
different reconstructions. If a real count rate dip in a cold trap were wrongly 
smoothed over in the reconstruction, then one would expect to see a radial 
variation in the mean residual moving out from the cold trap centre. The pat- 
tern of these spatially correlated residuals would depend upon the size of the 
cold trap, the depth of the count rate dip, and the FWHM of the instrumental 
response. 



3. 1 The radial dependence of the mean stacked residuals 



In order to try and discriminate between the two types of reconstruction, it 
is desirable to stack together the residuals around all sufficiently large (at 
least 3 contiguous pixels) cold traps. From the right hand panels in Fig.l, it is 
clear that the cold traps have a variety of sizes and shapes, and that they are 
often sufficiently close to one another that the anticipated rings of non-zero 
residuals around them would overlap. In the face of these complications, the 
easiest way to understand the significance of the results is to use a Monte Carlo 
method, whereby mock data sets are created and then analysed in exactly the 
same way as the real data itself. To this end, 100 different noisy mock Lunar 
Prospector time series data sets were created, taking the reconstructions in 
Fig.l as the input images, /(x). The different mocks varied only as a result 
of the random numbers chosen to add the noise. Both coupled and decoupled 
reconstructions were performed, and the stacked radial profiles of the reduced 
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residuals, r = R/cr with a coming from counting statistics, around cold trap 
centres were calculated. Pixels positioned near to more than one permanently 
shaded region were only counted once, at the radius appropriate for the nearest 
cold trap. 

The left hand panel in Fig. 2 shows the radial dependence of the reduced resid- 
uals for both types of reconstruction of the mocks created from the decoupled 
reconstructions of the real data. Error bars represent the la scatter between 
the results obtained from the individual mock data sets. The decoupled re- 
constructions show no radial trend in the mean reduced residual around the 
cold traps, which contrasts with the results from the coupled reconstructions. 
These show a negative residual for distances less than ~ 25 km, swapping 
over to a positive one at larger radii, signifying that I is too large in the cold 
traps. While the global \ 2 ma Y be acceptable, the coupled reconstruction pro- 
duces spatially correlated residuals near to cold traps. To remove these would 
require either a vastly increased number of pixons, which would produce a 
hugely detailed and improbable reconstruction, or the imposition of a shadow 
map. In short, these mock Lunar Prospector epithermal neutron data sets 
contain sufficient information to distinguish between the coupled and decou- 
pled reconstructions, i.e. the real data should tell us if the excess hydrogen is 
concentrated into the cold traps or more diffusely distributed. 

The right hand panel in Fig. 2 shows how the mean reduced residuals vary 
when the pixon algorithm is applied to mock data sets constructed from the 
input image where the count rate dips are not localised to the crater, i.e. the 
coupled reconstruction of the real data. It is apparent that there is a small 
bias, whereby the decoupled reconstruction still tries to put a slight count 
rate dip into the cold traps where it is not actually present. Also, the coupled 
reconstruction tends to oversmooth the small dips that did exist. However, 
these systematic biases resulting from the method are negligible relative to 
the size of the deviations found in the case of the decoupled input image. 

Fig. 3 shows the corresponding information for the coupled and decoupled re- 
constructions of the real data, with the error bars now representing the error 
on the mean reduced residual. A very similar pattern in the residuals is seen, 
for all but the central bin, as was the case for the mock reconstructions made 
using the decoupled input image. There are relatively few pixels, and hence 
observations, contributing to the central bin, which nevertheless still lies in 
the tail of the distribution of the mock results. Just as for the mock recon- 
structions, including the independently determined shadow map removes the 
correlated residuals out to 60 km from the cold trap centres. Thus, the Lunar 
Prospector data are not adequately described by the coupled reconstruction, 
and they do require that the epithermal count rate dips are concentrated into 
the permanently shaded regions. 
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3.2 Shadowiness 



The results in the previous subsection imply that the Lunar Prospector ep- 
ithermal neutron data themselves pick out the pixels that should be cold 
traps. In this section, the same conclusion is reached in a more quantitative 
way without resorting to stacking together the disparate cold traps. Consider 
a shadowiness parameter S defined via 



S = r * F, 



(13) 



where the filter function F, shown by the green curve in Fig. 2, is given by 



F(d) = 



0.035 tanh ( (d/k g~ 18 ) - 0.02 d < 60 km; 
d > 60 km. 



(14) 



F(d) is merely a function that represents the radial dependence of the spa- 
tial correlations in the reduced residuals when a coupled reconstruction is 
performed on a data set derived from an input image where the hydrogen is 
concentrated into shaded areas. For spatially uncorrelated residuals, the dis- 
tribution of S should have a median of zero. A significantly positive median 
value would imply that the reconstruction is systematically overestimating the 
count rate in that set of pixels. A map of the shadowiness parameter is shown 
in Fig. 4 for the coupled reconstruction of the real north pole data, where the 
shadow map prior is relatively well known. The locations of the independently 
determined shadow pixels are superimposed in white. Areas of high shadowi- 
ness, shown in blue, tend to coincide with the permanently shaded regions. 

Fig. 5 shows the cumulative distributions of shadowiness for a random reduced 
residual field, and for the shadow pixels for a number of different coupled 
reconstructions. The points show the mean distribution recovered from 100 
mocks made using either the decoupled (circles) or coupled (crosses) input 
north pole image. Error bars represent the standard deviation between the 
mock reconstructions. The long-dashed line is the distribution for the real 
north pole shadow pixels. If the shadowiness parameter and the shadow map 
were truly uncorrelated, then this cumulative distribution of S would be in- 
distinguishable from that resulting from the randomly distributed residuals. 
This is not the case, implying that the Lunar Prospector data have picked out 
as special, in a statistical sense, the pixels that are independently known to 
be permanently shaded. While not as convincing, perhaps as a result of the 
less well determined shadow map, the south pole results have also been added 
with a short-dashed line. 
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3.3 Uncertainties in the reconstructions 



While the stacked reduced residuals from all north and south pole shaded 
regions contain enough information to determine whether or not the count 
rate dips are typically concentrated into the cold traps, it is also interesting 
to know what is the uncertainty on the mean count rate within any individual 
cold trap. Permanently shaded regions that received more observations, either 
as a result of being large or near to the poles, will have better constrained 
mean epithermal neutron count rates. One would really like to know, given 
the measured cold trap count rates, what is the probability distribution of the 
true count rates. The mock catalogues allow two slightly different, yet closely 
related, questions to be answered, namely: 

1. if the hydrogen is as concentrated as the decoupled reconstruction of the 
real data suggests, then what are the bias and scatter in the decoupled recon- 
struction values, and 

2. if the hydrogen is as diffuse as the coupled reconstruction of the real data 
suggests, then what is the probability that a decoupled reconstruction will 
yield a cold trap count rate at least as low as did that from the decoupled 
reconstruction of the real data. 

Table .1 contains this information, calculated from the appropriate set of 100 
mocks for a few notable craters. From these data, it appears as if the north 
pole contains permanently shaded areas with deeper dips, i.e. more convincing 
hydrogen excesses, than exist in the south. An epithermal cou nt rate of ~ 13.5 



per s econd corresponds to 1 wt% water-equivalent hydrogen ([Lawrence et al. 
2006h . 



3.4 Discussion 



The results presented in this section demonstrate that: 

1. the Lunar Prospector epithermal neutron data rule out the hypothesis that 
the hydrogen concentration is smoothly distributed in the polar regions, and 

2. the set of perm anently shadowed pixels defined by the shadow maps of 



Elphic et al.l (120071 ) is statistically picked out as special. 



The first of these statements follows from the fact that spatially correlated 
residuals result if one assumes that the hydrogen distribution is smoothly dis- 
tributed throughout the polar regions. These statistically unacceptable residu- 
als are evident in both Fig. 3 and Fig. 5. Over most of the reconstructed image, 
the residuals are statistically acceptable and spatially uncorrelated. However, 
these figures show the problem because they focus on the residuals in the 
vicinity of a particular set of pixels, namely those that are deemed to har- 
bour regions of permanent shade. The information content latent within the 
data is sufficient to reject the hypothesis that the hydrogen is smoothly dis- 
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tributed, but one needs to know where to look to discover this message. The 
shadow map merely acts to focus attention into the regions where the useful 
information resides. If a random set of pixels were chosen as centres around 
which to investigate the residuals, then these correlations would not be seen. 
In fact, mapping the south pole shadow map pixel locations onto the north 
pole data and looking at the residuals around this different set of pixels shows 
no problem at all. The shadowiness distribution is completely consistent with 
that expected if the residuals were really random and spatially uncorrelated. 
Similarly, the residuals around a set of pixels in north pole craters supposedly 
without areas of permanent shade also show no deviations from the results that 
true randomly distributed residuals would produce. These observations illus- 
trate the second statement above. The areas that are independently thought 
to harbour permanent shade lie preferentially in regions where the smoothly- 
distributed hydrogen reconstruction produces spatially correlated residuals. 
The simplest, and only plausible, way to remove these local correlations in the 
residuals is to concentrate the hydrogen into the cold trap pixels. 



4 Conclusions 



This study shows that the Lunar Prospector data alone require the polar 
hydrogen excess to be concentrated into the permanently shaded cold traps, 
rather than being more diffusely distributed. In some craters the concentration 
of hydrogen is sufficiently high that it corresponds to ~ 1 wt% water. This 
combination of localisation into the permanently shaded polar craters and 
consequent relatively high concentration is consistent with the hypothesis that 
the hydrogen is present in the form of water ice. However, these concentrations 
are not sufficiently high to preclude the hy drogen being in the fo rm of solar 
wind protons implanted into regolith grains (IFeldman et all 120011 ) . While the 
Lunar Prospector data are not able to determine the intracrater distribution of 
hydrogen, radar results suggest that it is unlikely to be concentrated into small 
patches of high purity water ice, at least not in the regions of the cold traps 
accessible to Earth-based detectors. Assuming that the excess hydrogen is in 
the form of water ice that reac hes 2 m into the regoli th, the depth expected to 
be gardened in 2 billion years (IFeldman et all 1200 ll ) , this suggests that there 
are a few times 10 11 kg of water ice within permanently shaded regions within 
10° of the lunar poles. 
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Mean epithermal neutron count rates and their uncertainties for a few 

notable locations 



Crater 


Location 


d 


c 


dmock ± If 


P(< d\c) 


Cabaeus 


84.5°S,322°E 


14.8 


18.1 


14.1 ± 1.1 


0.04 


de Gerlache 


88.5°S,273°E 


16.8 


18.5 


17.5 ±0.8 


0.22 


Faustini 


87.3°S,77°E 


17.2 


18.4 


17.0 ±0.7 


0.22 


Shackleton 


89.7°S,110°E 


16.7 


18.8 


18.8 ±0.5 


< 0.01 


Shoemaker 


88.1°S,45°E 


17.8 


18.7 


18.4 ±0.4 


0.01 


Unnamed 


89.2°N,36°E 


10.5 


18.6 


10.0 ±1.3 


< 0.01 


Unnamed 


89.6°N,45°E 


13.8 


18.7 


12.8 ± 1.1 


0.01 



Table .1 

The decoupled, d, and coupled, c, reconstruction mean cold trap epithermal neu- 
tron count rates are listed, in counts per second, along with the mean and scatter 
of the decoupled reconstructions of the decoupled mocks, d m ock> and the probability 
of measuring a value at least as low as d, when doing a reconstruction of an image 
where the true value is c. 



14 



data 



coupled 



< m i ' I I | . .1 I J 1 I . I |i I I Ml I L I | I I I I 1 - I . | I I I I II 



100 ■ 



-100 



100 

I 

-100 E 



r l f I i ll 1 1 



i 1 1 1 1 1 1 1 1 H n 1 1 1 1 l ln m 1 1 1 1 1 ' i n 1 1 1 1 1 L 



.'.J l-f-l-^J.'i,-- i H j' i J ■ I- •■' 



"-SP 



I ■ .1 J I I I I ■ ■ ■ ■ ' ' I I I 



decoupled 

i i i I i i i ! — | i i i | | | | i; 



1 j i i i+i ■ j j I i i -i i j. i T i j 



so 



- IB 7 



a 

S- 



la | 

o 
u 



- 14 



u 



-100 100 -100 100 
x/km x/km 



-100 100 
x/km 



Fig. .1. Epithermal count rate maps for the north (top row) and south (bottom 
row) poles. From left to right, the columns represent the data, D, the coupled and 
decoupled reconstructions, I. The white circles represent a latitude of ±85° and 
the shading shows Clementine imaging. Also shown, with a black circle in the top 
left panel is the size of the Lunar Prospector's PSF. The spherical surface has been 
represented in these panels such that the distance from the image centre represents 
the arc length to the pole, and, following convention, 0° longitude points down/up 
for the north/south poles respectively. 
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Fig. .2. The variation of the mean reduced residual per observation, f, with the 
distance to the centre of the nearest permanently shaded crater, d. The results 
are stacked for all permanently shaded areas in both the north and south polar 
regions containing at least 3 contiguous shadow pixels. Coupled (black) and decou- 
pled (blue) reconstructions are shown in both panels. The left hand panel shows 
reconstructions from mock data made with the decoupled reconstruction of the real 
data as the input image. Error bars represent the la scatter between the individual 
mock results. The green curve represents the function F(d) given in Eq. 14, which is 
a fit to the results from the coupled reconstructions. The right hand panel shows the 
results for mocks made using the coupled reconstructions of the real data, shown in 
the central column of Fig.l, as the input image for the 100 mock catalogues. 
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Fig. .3. As for Fig. 2, using the stacked reduced residuals around cold traps for 
reconstructions of the real data. Error bars represent the errors on the mean reduced 
residuals. 
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Fig. .4. The shadowiness map for the coupled reconstruction of the north pole 
Lunar Prospector data. Cold traps are outlined in white, and these can be seen to 
lie preferentially in regions of high shadowiness according to the neutron data. 
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Fig. .5. Shadowiness distributions for various coupled north pole reconstructions. 
Points show mean results from 100 mock data sets made using either the coupled 
(crosses) or decoupled (circles) input image. Error bars represent the standard de- 
viation between the mock distributions. The distributions of the mock results are 
accurately represented by Gaussians. Results from the north and south pole recon- 
structions of the real data are shown by long- and short-dashed lines respectively. 
A solid line shows what the shadowiness distribution looks like for an uncorrelated 
random Gaussian residual field with a = 1. 
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